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Abstract 

The dynamical evolution of the refractive index of the tensor modes of the geometry pro¬ 
duces a specihc class of power spectra characterized by a blue (i.e. slightly increasing) slope 
which is directly determined by the competition of the slow-roll parameter and of the rate 
of variation of the refractive index. Throughout the conventional stages of the inflationary 
and post-inflationary evolution, the microwave background anisotropies measurements, the 
pulsar timing limits and the big-bang nucleosythesis constraints set stringent bounds on the 
refractive index and on its rate of variation. Within the physically allowed region of the 
parameter space the cosmic background of relic gravitons leads to a potentially large signal 
for the ground based detectors (in their advanced version) and for the proposed space-borne 
interferometers. Conversely, the lack of direct detection of the signal will set a qualitatively 
new bound on the dynamical variation of the refractive index. 


^Electronic address: massimo.giovannini@cern.ch 



1 Introduction 


It has been speculated long ago that gravitational waves might acquire an effective refractive 
index when they evolve in curved space-times [T]. As electromagnetic waves develop a refrac¬ 
tive index when they travel in globally neutral (but intrinsically charged) media, a similar 
possibility can also be envisaged in the case of linearized gravity. In this investigation it is 
suggested that the consistent variation of the refractive index throughout the conventional 
stages of the cosmological evolution leads to the production of a stochastic background of 
relic gravitons with blue spectral slopes. 

Relic gravitons are known to be produced in the early Universe thanks to the pumping 
action of the gravitational field [2]. This phenomenon occurs in a variety of different scenarios 
and, in particular, in the case of conventional inflationary models (see e.g. [a n El El d El 
El ED] for an incomplete but potentially interesting list of time-ordered references). While 
inflationary models typically predict decreasing slopes, in the conventional lore blue spectral 
indices of the relic graviton backgrounds can arise when a long phase (dominated by stiff 
sources) takes place after inflation but prior to the dominance of radiation [5]. Other less 
conventional possibilities include gravity theories which are not of Einstein-Hilbert type (see 
e.g. last paper of Ref. [ 3 ]) and the violation of the dominant energy condition in the early 
Universe m- In this paper we are going to argue that blue spectral slopes may arise from 
a comparatively more mundane possibility, namely the temporal variation of the refractive 
index of the tensor modes while the evolution of the background geometry follows exactly 
the same patterns of the concordance paradigm. 

The stochastic backgrounds of relic gravitons are subjected to three complementary 
classes of constraints. The first class of direct limits stems from the temperature and polar¬ 
ization anisotropies of the cosmic microwave background [ElllalEl] and it is customarily 
expressed in terms of a bound on the tensor to scalar ratio rxikp) at a conventional pivot 
wavenumber kp where the large-scale power spectra are assigned. The pulsar timing mea¬ 
surements [lEl ES] impose instead an upper bound on the cosmic graviton background at a 
typical frequency roughly corresponding to the inverse of the observation time along which 
the pulsars timing has been monitored. Finally the big-bang nucleosynthesis limits mi 
set an indirect constraint on the extra-relativistic species (and, among others, on the relic 
gravitons) at the time when light nuclei have been firstly formed. 

In this paper we intend to compute the cosmic background of the relic gravitons induced 
by the consistent variation of the refractive index during the early stages of the evolution of 
the geometry. The appropriately constrained spectra shall be compared with the frequency 
window of the ground-based interferometers such as Ligo/Virgo [HI ES], Geo600 [20] and 
the recently proposed Kagra [ 21 ] (ideal prosecution of the TamaSOO experiment [ 22 ] )• There 
also exist daring projects of wide-band detectors in space like the Lisa interferometer [ 23 ] 
(in one of its different incarnations) or the Bbo/Decigo [ 23 ] projedj^ 

^ The acronyms appearing in this and in the previous sentences refer to the corresponding projects: 
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The variation of the refractive index along the different stages of the evolution of the 
background must be continuous and differentiable at least once. This basic requirement 
stems directly from the evolution equations of the tensor modes of the geometry. In the case 
of a conventional inflationary and post-inflationary evolution the scale factor approximately 
evolves a^ 


a^(r) 

am(r) 



T < -n, 


I3t + {(3 + l)ri 

Tl 


-Ti < r < Ta, 


[/3r -f /3r2 + 2(/3 -1- 
4ri[/3r2 + {13 + l)ri] 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 


where ti coincides with the end of the inflationary phase and ra coincides with the time 
of matter-radiation equality; note that /? —)■ 1 in the case of a pure de Sitter phase and 
f3 = 1 — (9(e) in the quasi-de Sitter case. The transition to the domination of the dark 
energy will be discussed later on since, in practice, it does not affect the slope and it has 
a mild effect on the amplitude of the spectrum. Equations (inj. ( |1.2[ ) and ( |1.3[ ) are all 
continuous with their Erst derivatives at the transition point^ 

Since relic gravitons are produced because of the pumping action of the background 
curvature (containing second derivatives of the scale factor), the continuity of Eqs. (1.1), 
(1.2) and ( 1.3[ ) at the transition points ensures that the evolution equations of the relic 
gravitons will not have singularities in ri and T 2 but, at most, jump discontinuities. To 
guarantee a continuous variation of the refractive index without imposing further conditions, 
we are led to the following parametrization: 


n(r) = ?T.ia“(r), a > 0, 


(1.4) 


where a measures, in practice, the rate of variation of the refractive index in units of the 
Hubble rate. Equation (1.4) implies that n,(r) is automatically continuous and differentiable 
in Tl and T 2 provided the scale factor shares the same properties at the transition points. 
The parametrization (1.4) is minimal insofar as it contains only two arbitrary parameters, 
namely ni and a. Furthermore the value of rii is not totally arbitrarjj^ since at the onset 


Lisa (Laser Interferometer Space Antenna), Bbo (Big Bang Observer), Decigo (Deci-hertz Interferometer 
Gravitational Wave Observatory) and Kagra (Kamioka Gravitational Wave Detector). 

^We are assuming here a conformally flat Friedmann-Robertson-Walker background metric = 
a 3 {T)rj^i, where a(r) is the scale factor, r denotes the conformal time coordinate and is the Minkowski 
metric. This is the simplest way of complying with the concordance scenario where the extrinsic curvature 
is always much larger than the intrinsic (spatial) curvature. 

^This means, more specifically, that for r = —ti we have that off—ri) = ar{—Ti) and a'(— ti) = a(.(— ti). 
Similarly at the second transition point ar(T2) = am,{T2) and a(,(T2) = a'^(T2). 

^For practical reasons the bounds on the amplitude of the spectral index can be more simply expressed in 
terms of rii = n{Ti) where coincides with beginning of the inflationary phase. Since the evolution during 
inflation is known ni can be easily related to Ui by a simple redshift factor. 
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of the inflationary phase the refractive index must be larger than (or equal to) 1 to avoid 
superluminal phase and group velocities. Note that when the background is ever expanding, 
the positivity of a guarantees that this condition is preserved throughout the evolution of 
the geometry. 

The plan of this paper is the following. In section the most relevant technical aspects 
of the analysis are derived. The power spectra and the spectral energy density of the relic 
gravitons are computed in section and in the framework of the conventional cosmological 
evolution. In section]^ the spectra of the relic gravitons are confronted with all the available 
constraints. In the hnal part of section the prospects for the wide-band detectors of 
gravitational waves are illustrated. Section contains our concluding remarks. Some useful 
but lengthy results are collected in appendices and 


2 Relic gravitons with refractive index 


2.1 Basic definitions 

The two polarizations of the gravitational wave are dehned as 

= {niimj - ijij,), qf (^) = (Aift + ijiiri,), 


(2.1) 


where ki = ki/\k\, = mj/|m| and q = qi/\^ are three mutually orthogonal directions and 

k is oriented along the direction of propagation of the wave. It follows directly from Eq. 


(|2.l|) that ^ = 2(5 aa' while the sum over the polarizations gives: 


^ ) (jj/i (^) ^rnn(^) Pmiik^Pnjik^ T Pijik^Pmnik') 


( 2 . 2 ) 


where Pij{k) = {5ij — kikj). Dehning the Fourier transform of r) as 


(2.3) 


the tensor power spectrum Vrik^r) determines the two-point function at equal times: 

27r^ 

{hij{k,T) hran{p,r)) = j^Vrik, t) Sijmn{k)S^^\k + p) , (2.4) 


where, Sijmn(k) = J2xeij\k) e(,))^(fc)/4. The analog of Eq. (2.4) for h[j{k,T) is given by: 


^ 271 ^ 

{Kj{k,T) h'^^{p,T)) = —QT{k,T)Sijmn{k)S^^\k+p), 


(2.5) 


where Qrik^r) is the corresponding power spectrum and where the prime denotes a deriva¬ 


tion with respect to the conformal time coordinate r. Note that Eq. (2.4) follows the 
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same conventions used when deriving the spectrum of curvature perturbations on comoving 
orthogonal hypersurfaces (customarily denoted by 71{x,t)) 

{'}Z{k,T)n{p,T)) = -^VTi{k,T)S^^\k + ^, ( 2 . 6 ) 

which is exactly the quantity employed to set the initial conditions for the evolution of the 
temperature and polarization anisotropies of the Cosmic Microwave Background [laiiallll]. 
We also remind for future convenience that, according to the standard convention, the scalar 
power spectrum is assigned as: 

Vn{k) = Ani^Y) ’ A:p = 0.002 Mpc-\ (2.7) 

where kp is called the pivot scale, is the scalar spectral index and Ar is the amplitude of 
the scalar power spectrum at the pivot scale. 


2.2 Power spectra and spectral energy density 


The equation obeyed by hij{x,T) follows from the second order action: 




1 

M, 


(fx j dr a?{T) 


3t hijd-j- hij 


n" 


(^) 


3kh'ijd^hi 


(2.8) 


which reduces to the conventional action [101 ES] in the limit n(r) —)• 1. Note that in Eq. 


(2.8) ip = a/SttG = 8n/Mp and Mp = 1.22 x lO^^GeV. From Eq. (2.8) the equations of 
motion for hi are: 


h'k + 2Hh:, - 








= 0 , 


(2.9) 


where 77 = (In a)' = aH and H is the conventional Hubble rate. In Eq. (2.9) the contribution 


of the (transverse and traceless) anisotropic stress has been neglected. At low frequencies and 
in the concordance paradigm the contribution to the anisotropic stress is due to the presence 
of (effectively massless) neutrinos [29]. At high frequencies the anisotropic stress induced by 
waterfall helds may lead to an enhancement of the spectral energy density (see third paper 
in Ref. uni). Both effects will be neglected in what follows for two independent reasons. We 
shall neglect neutrinos because they are known to suppress the energy density of the relic 
gravitons at intermediate frequencies but their numerical relevance is not strictly essential 
for the present considerations. The waterfall held, on the contrary, may lead to large effects 
which are, however, model dependent, insofar as they arise in a given and specihc class of 
inhationary scenarios. 

In the absence of anisotropic stress hij{x, r) can be quantized and the corresponding held 
operator is: 

y/2ip 


hij{x,T) = 


(27r)3/2 ^, 


d^k e\j\k) [Fk,x{r)dj:^e 


-ik-x 


+ F, 


t Ak‘X] 


k,xA)a'^^e 


( 2 . 10 ) 
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where F^xir) is the (complex) mode funct ion obeying Eq. (2.9) and the sum is performed 
over the two physical polarizations of Eq. (2.1); note that 4 ,a'] The 

same expansion of Eq. (2.10) can be obtained for the derivative of the amplitude 

x/2ip 




(Fk elj’{k) [Gk,x{T)dj:^e 


—ik’X 


+ G 


k,X 


:T)d\ 


fcA 


( 2 . 11 ) 


where, this time, Gk = Fj^. The power spectra introduced in Eqs. (2.4) and (2.5) become, 
in this specihc case: 


VT{k,T) 

4/2 4,3 

= V 

(2.12) 

QT{k,T) 

4/2 4,3 

- V 

(2.13) 


The mode functions Tfc(r) obey the following equation which is the Fourier space analog of 


Eq. (Q: 


F'' + 2-F' + —F, = 0, 

a n^[T) 


(2.14) 


that can also be written as 


fk + 


k^ _ a 
.in?{T) a J 

Following Ford and Parker iza (see also [261 EZl 
energy density of the relic gravitons can be written as 


fk = 0. (2.15) 

for complementary approaches) the 


1 


Pgw — 


8£|,a2 L 


drhijd-rhij + 


77,2 (r) 


dkhijdkh. 


IJ 


(2.16) 


Within the established notation^ the energy density per logarithmic interval of wavenumber 
becomes 


dpgw 

dink 




k^ 

n^(r) 


Vrik^r) + Qrik.r) 


-)■ 


k^ 


4£p a^(r) n^(r) 


VT{k,T), 


(2.17) 


where the hnal result holds when the modes are inside the Hubble radius since, in this 
case, k’^Vpik^T)/rPir) —)■ QrikjT). In the opposite limit we have instead that Qrik^T) —)■ 
'H‘^VT{k,T). When discussing the graviton spectra over various orders of magnitude in 
frequency it is more practical to deal with the spectral energy density of the relic gravitons 
in critical units per logarithmic interval of wavenumber: 

1 dPgw 


(^5 


Peru = 3i7 /£ 


P- 


(2.18) 


Perit d In k 

The energy density of the relic gravitons per logarithmic interval of comoving wavenumber 
(or logarithmic interval of comoving frequency) introduced in Eqs. (2.17) and (2.18) will be 
occasionally called spectral energy density of the cosmic graviton background. 


®We take the opportunity for an elementary observation which is however rather crucial to avoid potential 
confusions: in this paper the natural logarithms will be denoted by “In” while the common logarithms will 
be denoted by “log”. 
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2.3 Practical time parametrizations 

We conclude this section with few remarks involving the time parametrizations. As we 
saw the evolution of the mode functions can be perfectly well discussed in the conformal 
time parametrization. However, for an explicit solution of the equations, it is convenient to 
use the 77 -time parametrization. Indeed, the action can be expressed in a simpler form by 
introducing a different time coordinate dehned by dr = n{r])dri. In this case the action of 


Eq. (|2.8|) can be expressed as: 

d^x / drjb'^^r]) 


S = 


1 


8il 


drjhijdrqhij 




6 ( 77 ) = 


aid) 

n{r]) 


(2.19) 


The mode expansion is analog to Eq. (2.10) and it is given by: 


>)) = E / ( 2 . 20 ) 


where, however, the evolution equation obeyed by Fk{jj) differs from Eq. (2.14) and it is 
given by 

^ 2 /^\^ 2 _ 

5772 + bKdrjJ drj ^ ^ ■ 

The evolution of the mode function rescaled through 6 ( 77 ) will then read 


( 2 . 21 ) 


ay* 

9772 


+ 


bKor]^/] 


fk = 0 , fk = Kv)Fk{,v)- 


( 2 . 22 ) 


The parametrization of Eq. (1.4) implies that power-law behaviours in the r-parametrization 


translate into power-laws in the 77 -parametrization. This is always true except for the case 
when the relation between 77 and r is logarithmic. This happens, for instance, when a = 1 
and the scale factor evolves during the radiation-dominated phase (i.e. Eq. (1.2)). The 


same thing happens when a = 1/2 and the scale factor is the one of dusty matter (as in 


Eq. (1.3)). Recalling Eq. (2.19) for the dehnition of 6 ( 77 ), if a = 1 we have that 6 ( 77 ) cx 77 
during the radiation epoch; similarly when a = 1/2 we also have that 6 ( 77 ) cx 77 during the 
matter phase. In these two cases Eq. (2.22) has a plane-wave solution in 77 . Even if the 


cases Qc = 1 and 0 = 1/2 must be separately treated, the results do not have a prominent 
physical meaning since they belong to a region of the parameter space which is anyway 
phenomenologically excluded. We shall therefore proceed in the discussion by assuming, for 
the sake of conciseness, that 07 ^! and a 7 ^ 1 / 2 . 


3 Power spectra in the different phases 


The evolution of the mode functions of Eqs. (2.14) and (|2.15[) must be solved by taking into 


account the evolution of the refractive index (see Eq. (1.4)) in each of the different stages 
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defined, respectively, by Eqs. (1.1), (1.2) and (1.3). At the practical level the strategy is to 


pass from the r-parametrization to the r^-parametrization and then transform back the ob¬ 
tained result in the conformal time coordinate. Since this procedure is algebraically lengthy 
but completely straightforward, we shall simply present the hnal result for the correctly 
normalized mode function and avoid pedantic details. We hnally mention that it is useful 
to introduce, in some of the forthcoming equations, the obvious notation 




k 


niatnf{ry 


Ur[T) = 


k 

nia“(r)’ 


uJmir) = 


k 

nia“(r)’ 


(3.1) 


where the scale factors in the different epochs are parametrized as in Eq. (1.1)-(1.3) 


3.1 Power spectrum during inflation 


When the scale factor and the refractive index are given, respectively, by Eqs. (1.1) and 


(1.4), the normalized solution of Eq. (2.15) is given by: 

V, 


fkir) = 






/i = 


2(1 — e )11 q ;/ 3 | 


(3.2) 


where H^^ygyr)] is the Hankel function of the hrst kind [3Ul 151] with argument 5 'j('r) and 
index p. Equation (3.2) has been derived in the case where the slow-roll parameter e = 
is constant in time. In this case it turns out that = 1/(1 — e) since aH = 
—1/[(1 — e)r]. In Eq. (3.2) the normalization \Vi\ = ^tt/ (2|1 -f a(3\) guarantees that, up 


to an irrelevant phase, Eq. (3.2) coincides with a plane wave in the large argument limit of 


Hankel functions. Since ojinfir) depends on r it is practical to introduce a single argument 

9iij) a^ 

air) = 

* \l + al3\ \l + ay\ni V ti) 


Inserting Eq. (3.2) into Eq. (2.12) we obtain, after some algebra, the explicit expression 


of the inflationary power spectrum: 


Vrik, r) = 8 


f Hi 


\kTi\ 


\MpJ 7r|l -f af3\ 


r \ 1 + 2/3 

nJ 




(3.4) 


where /i can also be expressed as p = (3 — e)/[2(l — e-f a)]. Using Eq. (3.3) and considering 


the modes that are larger than the Hubble radius, Eq. (|3.4|) becomes: 

ji^: 


Pr(i.fe„„) = C++(!;)/£-)’ 


^Since a > 0 and /? = 1/(1 — e) the absolute value is pleonastic. During the radiation and matter epochs 
the analog factors are not necessarily positive definite. To keep a homogeneous notation the absolute values 
have been always included even when not mandatory. 










































o6—nT(e,a) 

C(e,a) = — 


TT^ 



1 + 


a 

1 -e 


2—nx(e,Q) 


? 7 ,'r(e, a) = 3 — 


3-e 


(1 — e + a) 


(3.5) 


Denoting with the present valne of the critical fraction of radiative species (in the 
concordance paradigm photons and nentrinos) and with Ar the amplitnde of the scalar 
power spectrnm at the pivot scale (see Eq. (2.7)) the valne of kmax can be expressed, for 


instance, in Mpc ^ nnits: 


/ 


VMpc 


-1 


= 2.247 X 10 


23 


Hr 

Hi 


7 - 1/2 


0.01 


1/4 


2.41 X 10-9/ U.15 X lO-y 


1/4 


(3.6) 


In Eq. (3.6) 7 acconnts for the possibility of a delayed reheating terminating at an Hnbble 


scale Hr smaller than the Hnbble rate dnring inflation. 

In what follows, as already mentioned in the introdnction, we shall rather stick to the 
conventional case where the reheating is sndden and 7 = 1/2 (or Hi = Hr since the end of 
the inflationary phase coincides with the beginning of the radiation epoch). In more general 
terms, however. Hr can be as low as 10 “'^^Mp (bnt not smaller) corresponding to a reheating 
scale occnrring just prior to the formation of the light nuclei. 


In the limit a —0 and ni —)■ 1, Eq. (3.5) leads to the standard result, namely: 


lim VAKt) 

a-s>0, ni-5-1 TT \Mp 


k 


kr 


-2e 


(3.7) 


which implies Hi/M r = tt rp Ar/IQ with rp = 16e and ut = —rp/S. This kind of consis¬ 
tency relation (stipulating that the tensor scalar ratio exactly equals 16e) will not be valid 
anymore in the present context and the specific form of the tensor to scalar ratio will be 
used in section]^ to constrain the possible values of a. 

We observe that whenever a > 0 the spectral index can increase since it can be naively 
larger than the slow-roll parameter. More specihcally expanding np(e, a) in the limit e < 1 
we will have that 

(a — 2 ) ^ 


3a 

nr = -1" V-wrC 

1 -|- a (1-1- ay 


0{e^ 


(3.8) 


Two possible situations can be envisaged. If a > 1 the spectral slope in always violet 


q.e. 


sharply increasing); in the limiting case a S> 1 we have, according to Eq. (3.8), that np —)• 3. 
This possibility is strongly constrained by backreaction effects as we shall specifically see in 
section]^ If 0 < a < 1 the spectra are blue (i.e. slightly increasing) provided a > 2e/(3-|-5e); 
in the opposite case (i.e. a < 2e/(3 -|- 5e)) the conventional limit is recovered and np — —2e. 

The physical region of the parameters corresponds to the situation where at the onset of 
inflation the refractive index is larger than (or equal to ) 1. If this is the case a super luminal 
phase velocity is avoided. Indeed, denoting with Tj the initial time of the evolution we shall 
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have that0 


-aNt 

= j =nie *, 


Hi > 1 . 


(3.9) 


If n* = 1 we shall have that ni —)■ exp {aNt) where Nt is the total number of inflationary 
efolds. As we shall see the detailed discussion of section implies that rij must indeed be 
OiX) even if not strictly equal to 1. 


3.2 Power spectrum during the radiation epoch 


Following the conventions established in Eqs. (3.2) and (3.3) the expression of the mode 


function during the radiation epoch shall be written as: 
T)r 


fk{r) = , '' Jur{T)y{T,Ti) [c+{k,Ti) H^^'>[gr{T)] +C-{k,Ti)H^^'>[gr{T)]\, (3.10) 


where 


y{r,ri) = 


_ (/5 + 1)_ 


P = 


2 11 


iTPr I = 


TT 


a\ 


211 


a\ 


(3.11) 


In terms of ?/(r, ri) the argument of the Hankel functions gri^) is dehned as 

k 


9r{r) = - -|j/(r,ri) = 

1 — a ni 1 


a\ 


"r (/d + 1) 1^-" 

T H-- Ti 




(3.12) 


The coefficients c±(fc, ri) are complex and they obey |c+(/c, ri)p — |c_(/c, ri)p = 1. The exact 
expression of the two mixing coefficients is reported in Eqs. (A.l) and (A.2) of appendix [A| 
and can be determined by matching continuously the inflationary mode function (i.e. Eq. 


(3.2)) with the one of Eq. (3.10) in r = —ri; in formulae the following pair of conditions 
must be imposed: 


fri-ri) = fr\-n), 


df, 


(inf) 


dr 


df, 


(rad) 


dr 


(3.13) 


The requirements of Eq. (3.13) follow directly from the continuity of the scale factors and 


of the extrinsic curvature. We remind that the continuity of the scale factor guarantees, in 
the present approach, the continuity of the refractive index. The continuity of the extrinsic 
curvature (related to the conformal time derivative of the scale factor) guarantees that 
a"/a = 'h? + l-i' will have, at most, jump discontinuities. The exact expression of the mixing 
coefficients determined in the present situation reproduces the conventional results when 
a —?• 0 (see Eq. (A. 8) of appendix [A]) . 


®We recall that the parametrization of the scale factors given in Eqs. (1.1)-(1.3) stipulates that oi = 


a(—Ti) = 1. There are some who prefer to set Oq = 1 (where Oq is the present value of the scale factor) but 
this is not the convention adopted in the present paper 


10 





























The exact results of Eqs. (A.6) and ( |A.7 ) can be expanded in powers of ~ 1: 

c+{k,Ti) - C-{k,Ti) = 0{g^) + 0{g^), 


c+(/c,ri)+ c_(/c,ri) = 2 c_(/c, n) = 

gi g^r 


i£{a,e) + 0{g^) + C>((7j) 


(3.14) 


where g^. = gr{—'Ti) and g^ = gi{—Ti) (see also Eq. (A.3) of appendix [A|) . For practical 
reasons, the following combination 


£{ci, e) = 


|1 — al 


2^+^r(/i)r(p) _ 

81/5(1 — a) 7 r ^ |l + a 5 l 


{5[1 - 2(1 - a)p - 2pa] + 1 - 2/i}, 


(3.15) 


has been introduced in Eq. (3.14); note that 5(a,e) only depends on a and e since all the 
other auxiliary variables (i.e. p and p) are independent functions of a and e. The result of 
Eq. (3.14) can be made more explicit by using the expressions of g^ and gp to lowest order 
in xi = fcri, the approximate expression of |c_(xi)p is given by 

|c_(a;i)p = T^(q!, e) |1 + (3.16) 


Since Eq. (3.14) we have that c+(fc, ri) ~ c_{k, Ti), for r > —Ti the radiation power spectrum 
is therefore given by: 

Ak^fp 


Vrik^r) = 


~y{T,Ti)\c-{k,Ti)\ Jp[gr{.r)], 


(3.17) 


7r|l — a|a^(r) 

where Jp{gr) = [H^^\gr) + The argument of Jp[fi'r('r)] is gr (not g^) so that deep 

in the radiation epoch the power spectrum can be obtained in the limit gr S> 1. In this case, 
using the standard limits of the Bessel functions, the power spectrum becomes: 

2 


Vrik^T) = 


= 64/hfi Yfai 


71 \Mf 


u/r) |fcrih |c_(fc,ri)hcos^ [fi'r(n)]; 


(3.18) 


where the large argument limit of Jpigri^)] has been used. In Eq. (3.18) we can replace 
cos^ [gri^)] —)■ 1/2 as it is customary in this kind of analyses. Thus, from Eq. (3.18) we can 
also deduce the energy density and express it in critical units: 

ip, 




3.3 Power spectrum during the matter epoch 


(3.19) 


In the matter-dominated epoch, using Eq. (1.3) into Eq. ( |1.4[ ), the normalized solution of 
Eq. (2.15) for r > r 2 is given by: 

fk{r) = 


\J 2 iOJm{P ) 


\/wm(r)z(r,ri,r2)| d+(fc, ri, T 2 ) + d-{k, n, T2)Hl^\gm{T)] 


(3.20) 
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where 


/ N f 2(/3 + 1) \ 

Ti, Ta) = (^r + Ta +- - — -T^j , 




a = 


2 |1 - 2a|’ 


\T^m\ = 


TT 


211 - 2a\ 


(3.21) 


Notice finally that (7m('r) in Eq. (3.20) is dehned as 

UJmir) 


gmir) = 


|l-2a| 

k 


2(/5 + 1) 

r + Ta H- - — -Ti 

4ri[/?ra + (/? + l)ri] 


ni|l — 2a| 


(3^ 


r + Ta + 


2(/g + i) 

/S 


ri 


l-2o 


(3.22) 


As in the case of c±(fc,ri) also the mixing coefficients d±(fc,ri,ra) can be determined by 
continnous matching of the relevant mode fnnctions across ra. In this case we shall then 
impose 

^ p(rad) ^ p(mat) 

/rVj)=/fVj). 


Ot 


(3.23) 


The explicit form of (i±(fc,ri,ra) is reported in Eqs. (B.l) and (B.2) of appendix [B| together 
with a specihc discnssion of some relevant physical limits. 

As in the case of radiation the power spectrnm can be easily obtained in all the interesting 
regions. More specihcally, recalling Eqs. (B.l) and (B.2) we have that d+(/c,ri,ra) — 
d-{k,Ti,T 2 ) in the relevant physical limit. From Eq. (3.13) and (3.23) the power spectrnm 
dnring the matter phase is therefore given by: 


Vrik^r) = 


4kHl 


7r|l — 2Q;|a^(r) 


z(r,ri,ra) |d_(/c,ri,ra)| J^[gm{T)], 


(3.24) 


where Ja{gm) = [H^^\gm) + Hjj^\gm)]/‘2. Deep in the matter epoch the power spectrnm 
can be obtained in the limit gm ^ In this case, nsing the standard limits of the Bessel 
fnnctions, the power spectrnm becomes: 


Vrik^r) = 




TT \Mp 


nmir) \kTi\‘^ \d-{k,Ti,T2)\‘^ COS^ [9m{T)]. 


(3.25) 


In Eq. (3.25) we can replace cos^ [5'm('7‘)] —t 1/2 and we can also dednce the energy density 
and express it in critical nnits: 




(3.26) 


Using the same expansions discussed in the radiation case we can obtain, for instance, the 
leading order expression for \d-{k,Ti,T 2 )\‘^ for \kTi\ 1 and \kT 2 \ <C 1: 

/ ^ \ 2cx(p—cr) 

\d-{k,Ti,T2)\‘^ = M{a,e)l^-^j \kT2\‘^^'’~'^\ (3.27) 
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where A4{a,e) is given by: 


M(a,e) = 


22/^-9/?i+^r2(cr)r2(/i)[l + cr(l - 2a) + 2p(l - a)f 


— q/|2 |x _ 2a\^~'^^ |1 + a|3\^~‘^^^ 


X 


/3 + 1 


- (1 -a)p-p 


(1 + afi) 


n 2 


(3.28) 


and it is only function of a and e because the other parameters (i.e. p, p and a) are all 
independent functions of a and e and they have been dehned, respectively, in Eqs. (3.2), 


(3.11) and (3.21). 


An interesting limit of Eq. (3.27) is the one /3 —)■ 1 and a —)■ 0; in this limit p —>■ 


3/2, p —)■ 1/2 and a —)■ 3/2. In this case (setting also ni —)■ 1) we have that jd-f —>■ 
(9/64)|fcri|“^|A;r2|“^. This is the standard result for the mixing coefficients in the case of a 
transition from a pure de Sitter phase to the matter-dominated epoch passing through the 
conventional radiation dominances 

It is hnally possible to obtain a general expression encompassing the radiation and matter- 
dominated phases for the energy density of the relic gravitons in critical units. The expres¬ 
sions applies for modes inside the Hubble radius at the present time and it is given by: 

f Hq\ 2/3+a[l/6-(p-a)] / Hi \ 2-«[2(p-<t)+1]/2 / Hq\ «[2(p-<^)+l]/2 


X .^2(At-l-CT)-«gArta[2(^+cr)-a] 


p \ 4-2(/^+p) 


T {k,keq) - l + 6l(—j +Cl(^) 


P \2(p-cr) 


kn 


P \4(p-o-) 


(A:, kgq) , 


(3.29) 


eq^ 


o2a+3 

Mian) = - 

37r 


where hi and ci are numerical constants of 0(1); in Eq. (3.29) we have already expressed 


rii in terms of Uj, i.e. the initial value of the spectral index. The accurate value of hi and Ci 
can be also obtained numerically by computing the transfer function of the energy density of 
the relic gravitons introduced in (see e.g. pro]). Alternatively one can compute the transfer 
function for the power spectrum and then compute the energy density [1] . In both cases the 
idea is to integrate numerically the background and the mode functions across the matter- 


radiation transition. For the present ends what matters, however, is Eq. (3.29) in the limit 
k 3> keq. When a —?• 0 and rij —?• 1 and for k ^ k^q Eq. (3.29) reproduces the standard 
result 


VL 


SvrUeJ 


Hi 

Wp 


k 


kn 


-2e 


(3.30) 


®The limit of the exact expressions in this specific case (see Eqs. (B.6) and (B.7)) coincides with the limit 
of the general expression obtained above: this is a useful check of the whole algebraic consistency. 
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As in the standard case, when k 3> kmax we have that Vtg^ is exponentially snppressed as 
exp [—5k/kmax] where 5 is a numerical factor that can be estimated in a specihc model of 
smooth transition (see. e.g. [TO]). 


In Eq. (3.29) the ratio Hq/Mp is raised to an a-dependent power that disappears in the 
conventional limit of Eq. (3.30) (i.e. a —)■ 0). The explicit values of Hq/Mp and H^q/H q can 
be explicitly written as 


Mr 


= 1.228 X 10 


,-61 


ho 

or 


H. 


eq 


Hn 


= 10 


5.27 


fhln 


m) V 


hin 




VO.13647 V4.15 X 10-5 


- 3/2 


,(3.31) 


where ho is the indetermination on the present value of the Hubble rate, Hmo is the critical 
fraction of matter density and flijo is the critical fraction of radiation energy density. 

It is hnally useful to remark that the spectral slope of in the high-frequency branch 
(i.e. for k 3> k^q) is simply given by 4 — 2(/i -|- p) as it can be immediately verihed from the 
explicit expression of Eq. (3.29). Recalling Eqs. (3.2) and (3.10) the high-frequency slope 


can be written more explicitly: 


4 - 2[p(e, a) + p{a)] = 3 + 


+ r 


a 


a 


1 -|- cr 


-~2(a-e) + C>(e2) + C>(a2)^ (3.32) 


where the second equality follows by expanding the exact expression Erst in powers of e and 
then in powers of a. This limit captures an important corner of the parameter space (see 
the discussion of section]^. Equation (3.32) implies that in the pure de Sitter limit without 
variation of the spectral index (i.e. a = 0 and e = 0) Qgw is constant in frequency with 


amplitude given by Eq. (3.30) in the limit k ^ k, 


''eq- 


3.4 Typical frequencies 


We shall always use wavenumberTH expressed either in units of Mpc ^ or in units of Hz. 


The reason of this potential ambiguity is that the discussion mixes constraints arising over 
large length-scales (where the wavenumber are typically measured in Mpc“^) and other 
limits coming from comparatively much shorter scales (where the wavenumbers ate typically 
assigned in Hz) [TU]. For instance, in what follows we shall be dealing with the big-bang 
nucleosynthesis wavenumber kbbn 


kbbr,. = 1.47 X 10-^° 


9p 

10.757 


Tbbn \ 

MeV7 


hln 


RO 


4.15 X 10-5 


1/4 


Hz, 


(3.33) 


where Pp denotes the effective number of relativistic degrees of freedom entering the total 
energy density of the plasma and Tbbn is the big-bang nucleosythesis temperature determining 
the size of the Hubble radius at the corresponding epoch. The typical value of the frequency 

shall often measure comoving wavenumbers in Hz and refer to typical comoving frequencies. Note 
that, in natural units, k = Frequencies and wavenumbers are not exactly coincident even if it is useful, 
at a practical level, to measure wavenumbers in Hz. 
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corresponding to Eq. (refkbbn) is Phbn = kbbn/‘^'^ = 2.3 x 10 Similar observations can be 
made in all the other cases. For future convenience k^ax and keq can also be expressed in 
Hz: 


kr^ 


= 2.183 


~h) 


0.01 


'/v 


V2. 


41 X 10-9 


) u 


15 X 10 


-5 


h — 


9.69 X 10-1" { -1/2 

vo.i364y V4.15 X lo-sy 


Hz. 


GHz, (3.34) 
(3.35) 


The frequencies corresponding to the hducial values of the parameters given in Eqs. (3.34) 
and (3.35) are given, respectively, by Umax = 0.34 GHz and by u^q = 1.54 x lO-i" Hz. 


3.5 Secondary effects 

In the present analysis we neglected, for the sake of simplicity, a number of secondary effects 
that may interfere with the variation of the refractive index. For k < k^tn the power spectra 
and the energy density of the gravitons are suppressed due to the neutrino free streaming. 
The effective energy-momentum tensor acquires, to hrst-order in the amplitude of the plasma 
fluctuations, an anisotropic stress (see e. g. [29] and references therein). The overall effect 
of collisionless particles is a reduction of the spectral energy density of the relic graviton^^ 

The second effect leading to a further suppression of the energy density is the late dom¬ 
inance of the dark energy. The redshift of A-dominance is given by Gde/GMo- In principle 
there should be a break in the spectrum for the modes reentering the Hubble radius after ta- 
This tiny modihcation of the slope is practically irrelevant and it occurs anyway for k < keq- 
However, the adiabatic damping of the tensor mode function across the TA-boundary reduces 
the amplitude of the spectral energy density by a factor (Gmo/I^a)^ — 0.10. This hgure is 
comparable with the suppression due to the neutrino free streaming. These effects have been 
discussed in the past (see mm and references therein). Further effects leading to similar 
reductions of flg^ are related to the evolution of the relativistic species. 

The effects mentioned in the two previous paragraphs are secondary since they can be 
easily reabsorbed by the variation of one of the other unknown parameters of the cosmic 
graviton background. At the same time they become truly essential if the absolute nor¬ 
malization of the graviton spectrum is known (see last paper of Ref. [IHI)- In the present 
case the inclusion of these secondary effects is unimportant for the hnal conclusions, as we 
explicitly checked. 

Assuming that the only collisionless species in the thermal history of the Universe are the neutrinos, 
the amount of suppression can be parametrized by the function F{Ru) = 1 — 0.539i?i, -t where R^, 

is the fraction of neutrinos in the radiation plasma. In the case N,^ = 3, = 0.405 and the suppression 

of the spectral energy density is proportional to (0.405) = 0.645. This suppression will be effective for 
relatively small frequencies which are larger than keq and smaller than khtn- 
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4 Phenomenological considerations 


The limits on the variation of the refractive index over various scales will now be derived. 
There are four qualitatively different sets of bounds to be examined and they involve, re¬ 
spectively, (i) the backreaction constraints during inflation, (ii) the limits stemming from 
the tensor to scalar ratio obtained from the temperature and polarization anisotropies of 
the cosmic microwave background, (in) the bounds arising from the millisecond pulsar tim¬ 
ing measurements and finally (iv) the so-called big-bang nucleosynthesis constraints. At 
the end of the section the impact of the derived limits on the prospects for the wide-band 
interferometers shall be addressed. 


4.1 Limits from backreaction effects 


The considerations of section suggesting an upper limit on a can be made more concrete 
by computing the total energy density of the gravitational waves and by comparing it with 
the critical energy density during inflation. From Eqs. (2.17) and (2.18) the total energy 


density of the produced gravitons in critical units is given by: 


Pgw ifli 0^5 
pcrit 


1 rVn dk 

24Hfa^Ji/n T 


n2(r) 


Vrik^r) + Qxik.r) 


Pcrit 


SHfMj 
Stt 


(4.1) 


where the integration is extended from modes the exiting the horizon at the onset of inflation 


up to those reentering exactly at the onset of the radiation phase. Inserting Eqs. (3.2), (3.3) 


and (3.4) into Eq. (4.1) and performing the indicated integrals we can easily obtain the 


following result 
(a, a, e) 


Pgv. 


pcrit 


= iS(a,e 


(H, 


\Mp 

{3-2ii)/0 


2fia 


ai\(3-2/.)//3r 
a J 


1 + 


(3 - 2/i) 


1 + 


(3 — 2/i) / a 


5(a,e) = 


(5 - 2/i) nf \ai 
2^^r^(/i) |l + a/jp^-^ 

3 TT^ (3 — 2/i) nf 


-2a 


(5- 

2/P' 


2/i) nf 


-2a 


a J 


(4.2) 


The function 5(Q:,e) appearing in Eq. (4.2) only depends on a and e since (i = /9(e) and 


/i = /i(a, e); moreover, the dependence on the scale factor in Eq. (4.2) can be traded for the 
total number of inflationary efolds Nf 

It is not necessary to analyze the independent variation of Nti cn and e: the upper bound 
on a is anyway less constraining than the ones to be examined later on. In fact a cannot 


^^The cases /i = 3/2 and p = h(2 are singular: this simply means that the corresponding integrals must 
be separately computed and lead to a logarithmic contribution which is only present, strictly speaking in 
the case e —)■ 0 and a —)■ 0 (i.e. pure de Sitter evolution). 
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exceed 0.1 when the remaining parameters are hxed to their hducial values: from Eq. (4.2) 
with a = ttf = ai and rii = 0{1) we have that \pgwlPcrit\ < 1 provided 


a < 


In (Tre^- 


n) 


3iV, 


(4.3) 


where, as usual, An is the amplitude of the scalar power spectrum at the pivot scale and has 
been introduced in Eq. ( |2.7 ). The total number of efolds Nf appearing in Eq. (4.5) must be 
larger than (or equal to) N^ax- 


Nmax = 61.43 + - In 






4.15 X 10-5 


-ln(^)+^ln 


A 


-R 


2.41 X 10-9 


+ ^ In 
4 


0.01 


(4.4) 


which is the maximal number of efolds presently accessible to large-scale observation^^^j32] . 
In the case (e, Nt, An) = (0.01, 65, 2.41 x 10-9) Eq. (4.4) implies, for instance, a < 0.11 
when rij = 0{1). Similar results can be obtained from slightly different choices of parameters. 


Nt = 60, k- 0.002 Mpc“', e- 10"^ Nt = 65, k- 0.002 Mpc"', e= 10“^ 



Figure 1: The value of tt computed from Eq. (4.5) is illustrated in the (a, n^) plane. 


4.2 Limits from the tensor to scalar ratio 

The long wavelength gravitons induce direct temperature and polarization. Technically they 
can affect the TT power spectra (i.e. the temperature autocorrelations) the EE power 

practice N^aax is determined by redshifting the inflationary event horizon at the present time and by 
identifying the obtained results with the current value of the Hubble radius. 
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spectra (i.e. the polarization autocorrelations) and the TE power spectra (i.e. the cross¬ 
correlation between temperature and polarization). These power spectra can interfere with 
temperature and polarization power spectra induced by the scalar mode and this is why the 
upper bounds on the tensor to scalar ratio can be derived from the accurate determinations 
of the temperature anisotropies and polarization anisotropies [T^l [13] • Depending on the 
combined data sets the WMAP 5-year data provided bounds on in the framework of the 
concordance model with values ranging from tt < 0.58 to tt < 0.2. Similar bounds have 
been obtained from the WMAP 7-year data. The 9-year WMAP data release gave a limit 
Tt < 0.38 always in the light of the concordance model in the presence of tensor. In the last 
three years there have been more direct determinations of the i?-mode polarization of the 
cosmic microwave background. The hrst detection of a B-mode polarization, not caused by 
relic gravitons but coming from the lensing of the E-mode polarization, has been published 
by the South Pole Telescope [33]. The Bicep2 experiment [33] claimed the observation of a 
B-mode component with = 0.2tQ;o5 which turned out to be induced, at least predominatly, 
by a polarized foreground. The present Planck data imply < 0.1 na but the reported 
sensitivity to the B-mode polarization is rather poor. 

While tensor contribution to the cosmic microwave background observables can be cleanly 
ruled out (or ruled in) by direct observations of the B-mode polarization (as attempted by 
Bicep2 and by other previous experiments directly sensitive to polarization (see e.g. [33])) 
for the present ends what matters is not the specihc value of the bound but the generic order 
of magnitude that tt should not exceed at a conventional pivot scak^^ kp = 0.002 Mpc“^ 
From Eq. (3.5) the tensor to scalar ratio reads 


rT{kp,a,ni,e,Nt) = vr e C{e,a) ^2a[3 nTie,a)]Nt ^ 


k \nT{e,a) 


kn 


(4.6) 


To get a superhcial idea of the orders of magnitude involved we can hrst consider the case 
Nt = 65 and e = (P(10“^). In this case we have that < 0.1 provided 0 < a < 0.13 for 
Ui = 1. A slight increase of n* or of the total number of efolds strengthen the limit on a. 
For instance if n* = 1.5 we will have that 0 < a < 0.04 (for W = 65) and 0 < a < 0.02 
(for Nt = 70). A more detailed discussion that is summarized Figs. [^and|^ In Fig. [^all 
the parameters are hxed except rij and a. Along each of the curves logr-r is constant and 
the labels refer to the value of the common logarithm (i.e. to base 10) of tt computed from 
Eq. (4.5). For illustration the arrows indicates the curve logr-r = —1: the physical region, 
compatible with the current constraints, demands that \ogrT < —1. 

In Fig. (plot on the left) the total number of efolds is W = 60, while in the plot on the 
right Nt = 65. As the number of efolds increases the value of a is pushed towards 0.1. 

The same trend is observed in Fig. (plot on the left) where is illustrated in the 
{a, Nt) plane. As the total number of efolds increases beyond 65, a is driven towards 0. 
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The WMAP collaboration consistently chooses kp = 0.002 Mpc ^. The Bicep2 collaboration used kp = 


0.05 Mpc ^. The first data release of the Planck collaboration assigned the scalar power spectra of curvature 
perturbations V-ji at kp = 0.05Mpc“^ while the tensor to scalar ratio rx is assigned at kp = 0.002 Mpc“^. 
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n,^l, k= 0.002 Mpc-^ e= 10"^ n;=l, Nt=65, k- 0.002 Mpc“i 



Figure 2: The value of is illustrated in the (a, Nf) plane and in the (a, loge) plane. 

Whenever gets smaller than 1 the parameter space in the {a, ni) plane gets larger, 
depending on the value of n*. This region has been excluded since it would lead to a 
superluminal phase velocity which coincides, in this case, with the group velocity. In specihc 
situations where the group velocity does not coincide with the phase velocity, the regions 
0 < Uj < 1 might become phenomenologically viable. However in the present context we 
just want to focus on the most conservative situation. 


4.3 Limits from the pulsar timing bound 

The bounds on tt are derived in the hypothesis that the consistency relations between the 
tensor amplitude and the tensor spectral index are verihed. This is not necessarily true in 
the present context. Thus the bounds stemming from might be even less stringent than 
the ones we just analyzed. Ultimately this is not an important limitation since the most 
constraining bounds, in the present situation, come from higher wavenumbers (or higher 
frequencies). Indeed, the pulsar timing constraint demands 

^{kpulsar,ro) < 1.9 X 10"®, kpulsar ^ 10"^ Hz, (4.6) 


where kpui^ar roughly corresponds to the inverse of the observation time along which the 
pulsars timing has been monitored. The same strategy discussed in the previous subsection 
can now be applied to Eq. (4.6). In Figs. and we used exactly the same range of 
parameters already employed in Figs. [^and[^ 

In Figs. [^and|^we illustrate the common logarithm of computed from Eq. (3.29). 
Comparing Figs. with Fig. the values of a allowed by the pulsar bound are much 
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0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 


a 


a 


Figure 3: The value of the common logarithm Vtgy^ computed from Eq. 
in the {a, rii) plane and at the pulsar frequency. 


(3.29) is illustrated 


smaller than 0.1 for the same range of variation of Uj. The same conclusion follows from the 
comparison of Figs. |^and[^ The pulsar bound is systematically more constraining because 
when the refractive index is dynamical the energy density of the relic gravitons is increasing 
(rather than decreasing) for k ^ keq- In Figs. and the arrow indicates, approximately, 
the curve where the pulsar bound is saturated. 


4.4 Limits from the big-bang nucleosynthesis 

A conclusion compatible with the pulsar bound can be drawn in the case of the big-bang 
nucleosynthesis constraint (BBN in what follows) stipulating that the bound on the extra- 
relativistic species at the time of big-bang nucleosynthesis can be translated into a bound on 
the cosmic graviton backgrounds. This constraint is customarily expressed in terms of 
representing the contribution of supplementary neutrino species but the extra-relativistic 
species do not need to be fermionic. If the additional species are relic gravitons we have: 


hi 


^bbn 


^Gw{k,To)d\nk = 5.61 x 10 


^0^70 


2.47 X lO-'^ 


(4.7) 


where k^bn and kmax have been computed, respectively, in Eqs. (3.33) and (3.34); note that 
in Eq. (4.7) denotes the present critical fraction of energy density coming just from 


photons. The bounds on ANy range from ANy < 0.2 to ANy < 1. The bounds stemming 
from Eq. (4.7) can be easily inferred from Figs. and In both cases we illustrate the 
common logarithm of the left hand side of Eq. (|4.7[) computed in the case of a dynamical 
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a a 


Figure 4: The common logarithm of at the pulsar frequency is illustrated in the (a, N^) 
plane and in the (Q!,loge) plane. 


N, = 60, e= 10“^ 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 


a 


a 


Figure 5: The common logarithm of the big-bang nucleosynthesis bound computed from Eq. 


(4.7) is illustrated in the (a, rii) plane. 


refractive index from Eq. (3.29). By comparing Figs. [^and|^we see that the pulsar and the 


BBN bound are largely compatible. Conversely by comparing Figs. and the big-bang 
nucleosynthesis bound is always more constraining than the bounds stemming from rrp{kp). 
The same conclusion is reached if we compare the variation of a different set of parameters, 
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a a 


Figure 6: The common logarithm of the big-bang nucleosynthesis bound computed from Eq. 
(4.7) is illustrated in the (a, Nt) plane and in the (Q!,loge) plane. 


like in Figs. ii and 1^ Note hnally that in Figs. and the arrow has been used to 
guide the eye towards the curve where the big-bang nucleosynthesis bound is approximately 
saturated. 


4.5 Prospects for wide-band detectors 

The bounds examined in the previous subsections suggest that for the standard hducial 
values of {Nt-, e) the values of rii and a are constrained to be within the following window: 

1 < ui < 10, 0 < a < 0.07. (4.8) 

The sensitivity of a given pair of wide-band detectors to a stochastic background of relic 
gravitons depends upon the relative orientation of the instruments (see e.g. [SniEZ!) and a 
specihc analysis of the signal to noise ratio is beyond the scopes of this paper. While the 
advanced version of the wide-band interferometers is still matter of debate (and the published 
results do not seem conclusive at the moment) the frequency window of the detectors will 
always be between few Hz (where the seismic noise dominates) and 10 kHz (where the shot 
noise eventually dominates). The wideness of the band (important for the correlation among 
different instruments) is not as large as 10 kHz but much narrower. There are projects 
of wide-band detectors in space like the Lisa, the Bbo/Decigo. The common feature of 
these three projects is that they are all space-borne missions; the Lisa interferometer should 
operate between 10“^ and 0.1 Hz. Nominally the Decigo project will be instead sensitive to 
frequencies between 0.1 and 10 Hz. 
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Recalling the results obtained so far we can safely say that growing spectra can arise in 
the following range: 

(3+5e) 3Af, ’ ^ ’ 

where the lower bound of Eq. ( 4.9[ ) has been derived after Eq. (3.8) while the upper bound 
follows from Eq. (4.3). Since the the upper limit of Eq. (4.9) is larger than the one of Eq. 
(4.8) we can conclude that growing spectra arise in practice in the whole range of variation 
of a. Using the illustrative example discussed before we have that for kusa = 0(10“^) Hz 
(i.e. compatible with the Lisa window) we would have Hgw(mHz) ~ for a = 0.06, 

Nt = 65, e = 0.001 and rii = 1. Larger values of n* and Nt make the signal smaller. For the 
Ligo/Virgo frequencies and for the same parameters chosen in the Lisa case we would have 
instead Hg„(0.1kHz) ~ 10“^'®® where kLigo = 0.1 kHz. 

This trend is conhrmed by the results illustrated in Fig. where we report the common 



Figure 7: The energy density of the relic gravitons is illustrated in the {a, rii) plane in the 
Lisa window (plot on the left) and in the Ligo/Virgo window (plot on the right). 


logarithm of in the Lisa window (plot on the left) and in the Ligo/Virgo window (plot 
on the right). The dashed lines in both plots corresponds to the common logarithm of the 
big-bang nucleosynthesis bound in the AN^ = 1 case (i.e. log (5.61 x 10“®) = —5.25). The 
allowed region of the parameter space must be, in both plots of Fig. below the dashed 
line^^ While the noise power spectra of Lisa are still rather hypothetical, the advanced 
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^By comparing the two plots of Fig. [^the dashed line (representing the big-bang nucleosynthesis bound) 
looks closer to the actual value of the logflgu, in the Ligo/Virgo case than in the Lisa case. Indeed, if the 
integrand increases, the integral appearing in Eq. (4.7) can be approximated by the value of flgw at kmax- 
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version of terrestrial interferometers might get down to 10“^° in after an appropriate 
integration time. The shaded area in Fig. [^illustrates the allowed region in the parameter 
space where the relic gravitons are potentially detectable. As already mentioned we shall 
not dwell here on the detectability prospects in this paper. The value = 0(10“^°) has 
been mainly quoted to guide the eye and whether this sensitivity will be in fact reached by 
terrestrial of space-borne detectors is an entirely different issue. If this value will be indeed 
reached by cross-correlation of two instruments more detailed discussions could be necessary 
(see e.g. [36l EZ] and discussions therein). 

It is hnally appropriate to mention that the maximal signal due to the variation of the 
refractive index occurs in a frequency region between the MHz and the GHz. In this range 
of frequencies microwave cavities or even wave guides [SHlEHlllo] can be used as detectors 
of gravitational waves. 


5 Concluding remarks 


The continuous and differentiable evolution of the refractive index leads to power spectra 
and spectral energy densities of the relic gravitons that are slightly increasing as a function 
of the comoving wavenumber (or of the comoving frequency). While the rate of variation 
of the refractive index can be stringently bounded, the derived limits do not exclude the 
potential relevance of the cosmic graviton background for either ground based or space- 
borne interferometers aimed at a direct detection of gravitational waves. The spectral slopes 
are determined by the competition of the slow-roll parameter against a which measures the 
rate of variation of the refractive index in units of the Hubble rate. The phenomenologically 
allowed region implies that 1 < n, < 10 and 0 < a < 0.07, where n* denotes the value of the 
refractive index at the onset of the inflationary expansion. 

The sensitivity of wide-band detectors to the relic graviton backgrounds is customarily 
expressed in terms of the minimal detectable spectral energy density of the relic gravitons 
in critical units. The operating windows of the ground based and space-borne interferome¬ 
ters are complementary: the typical frequency range of space-borne interferometers extends 
between 0.1 mHz = lO^'^Hz and few Hz. Conversely the window of ground based detectors 


i6We recall, for the sake of precision, that the expression of the signal-to-noise ratio in the context of 
optimal processing required for the detection of stochastic backgrounds [MlEZl depends on an integral over 
the frequency band (between few Hz and 10 kHz in the case of the ground-based interferometers of Ligo- 
type). The numerator of the integrand contains while the denominator we have the sixth power of the 
frequency multiplied the noise power spectra of each of the two correlated detectors. The signal-to-noise 
ratio depends also on the relative orientation of the interferometers and on the total observation time which 
is crucial to increase the sensitivity. Naively, if the minimal detectable signal (by one detector) is h^Q.gyj, 
then the cross-correlation of two identical instruments might increase the sensitivity by a factor Ij^fKiTT 
where /S.v is the bandwidth and T, as already mentioned, is the observation time. Therefore if a single 
instrument detects ~ 10“® the correlation may detect /iqHgw — 10“^° provided ~ 100 Hz and 

T ~ (P(lyr) [TO] . 
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V 

log 

log [5/iHz] 

log [a/^/iHzJ 

10“^ Hz 

-8.15 

-32.56 

-16.28 

10-2 Hz 

-7.94 

-38.35 

-19.17 

1 Hz 

-7.74 

-44.14 

-22.07 

102 Hz 

-7.53 

-49.93 

-24.96 

10^ Hz 

-7.32 

-55.73 

-27.86 

10° Hz 

-7.11 

-61.52 

-30.76 

10® Hz 

-6.90 

-67.31 

-33.65 

10^° Hz 

-6.69 

-73.10 

-22.07 


Table 1: Typical values of the spectral energy density and of the strain amplitude in the op¬ 
erating windows of ground-based and space-borne interferometers. Consistently with the dis¬ 
cussion of section|^the parameters {ui, a, e, Nt) have been chosen to be (1, 0.06, 0.001, 65). 


extends between few Hz (where seismic noise dominates) and 10 kHz = 10^ Hz (where shot 
noise dominates). While the time scales for the realization of space-borne interferometers 
are still vague it is useful to illustrate our hndings by keeping ground based detectors and 
space-borne interferometers on equal footing. 

In Tab. the frequencies encompass the operating ranges of space-borne and ground- 
based detectors. In the hrst column we illustrate the common logarithm of flgw In the 
two remaining columns we illustrate the common logarithms of the strain amplitude Sh and 
of its square roolj^ In Tab. [T] the range of frequencies has been extended well beyond 
the window of ground-based interferometers. The last frequency, i.e. 10 GHz is even larger 
than Umax = ^maa:/(2vr) = C>(GHz) and we just included it to show that the spectral energy 
density is always well below the constraints previously discussed. 

If the correlation of (advanced) wide-band interferometers will eventually reach sensitiv¬ 
ities 0(10“^°) in Qgu), the present considerations might become relevant since the signal due 
to a dynamical refractive index with a = 0(0.06) and n* = 0(1) can even be 0(10“^'®^) for 
the phenomenologically allowed region of the parameter space and for a typical frequency 
0(0.1 kHz). Even if the noise power spectra and the specihc features of the hypothetical 
space-borne interferometers are still unclear, a potential sensitivity 0(10“®'°®) in Qg^ cannot 
be excluded for a typical frequency O(mHz). The lack of detection of a stochastic back¬ 
ground of relic gravitons either by ground based detectors or by space-borne interferometers 
may therefore provide a further (and potentially much more stringent) bound on the rate of 
variation of the refractive index. 


^^The sensitivity is often expressed by means of Sh or in terms of its square root[l0]. The precise relation 
between Sh and flgu, is given by Sh(v,TQ) = 7.981 x 10““^^ (lOOHz/i/)^ hQ^g^{iy,To) Hz“^. Note that Sh is 
measured in 1/Hz = sec. 
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A Mixing coefficients during the radiation phase 

The explicit expression of the mixing coefficients appearing in Eq. ( 3.10[ ) can be written as 
c+{k,Ti) = , (A.l) 

c_(A;,ri) = (A.2) 


where q = ^|1 — q;|/|1 + a(3\. In Eqs. (A.l) and (A.2) the following notation have been 
employed: 

9i{-ri) =gi= gr{-ri) =gr= A-S) 


ni|l + a(5\ ’ 


ni|l — a|/9 


Furthermore, in Eqs. (A.l) and (A.2) and are two auxiliary 


defined as: 


= H^^\g,) 

= H^\gd 


(1 -a)p+ - 

(1 + a(3)g, + - 


- g^{l - a)H^pl^{g^), 


expressions 

(A.4) 

(A.5) 


where, following the properties of the Hankel functions under complex conjugation, we shall 
have = A^p^*{g^) and In more explicit terms Eqs. (A.l) and 

(A.2) can be written as: 


c+(fc,ri) = 


|1 — a| 






4(1 — a) \ |1 + a/3| 

(1 - a)j, 5. 


C-{k,Ti) = 




4(1 — a) \ |1 + a(3\ 


a\ 




P 

/3 + 1 

~w 


(A.6) 


+ (1 — ol)p + 


p(l + Qf/S) 


/S 


(1 - a)9, apmafUaP - * Ef(a,) <1(9,)}. 


(A. 7) 


As it can be explicitly verified from Eqs. (A.6) and (A.7), |c+(/c,ri)p — |c_(/c,ri)p = 1. It 
is useful to mention that Eqs. (A.6) and (A.7) reproduce exactly the standard results in 
the limit a —?• 0, rii —)■ 1 and /3 —)■ 1. More specifically this limit refers to the situation 
where there is a transition from an exact de Sitter phase (i.e. e —)■ 0 and /r —?• 3/2) to a 
conventional radiation-dominated phase where p —)■ 1/2. In the standard limit we also have 
that Pj —)■ kri and g^ —)■ kri and, from the explicit form of Eqs. (A.6) and (A.7), 

i 


c+{xi) = ^^[2xi{xi +i) - l]e 


2ixi 


C-{xi) = 


2xr ^ ^ ' 2xY 

in this case the mixing coefficients depend on the single argument Xi = kxi. 


(A.8) 
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B Mixing coefficients during the matter phase 

The explicit expression of the mixing coefficients appearing in Eq. (|3.20[) can be written as 


d+{k,Ti,T2) = 


m 


4^2(1-2a) \ |1 


|1 -2a| 


a 


+ C_{k,Ti) 
d_{k,Ti,T2) = 


m 


4^2(1-2a) ^ |1 


|1 -2a| 


a 


c+(«:,Ti) (9.) 

(B.l) 

- 2^^?>(9».)6f (9, 


C+(k,Ti 


+ C_(k,Ti) 




(B.2) 


where gm and gr are dehned as: 


QmiT'i) 9m 


2kT2 


f—V 

ni|l — 2a| V/3r2/ 


1 + 


/3 + 1 /^'Tl 


T2 


-| 1—a 


QriT'i^ 9r 


11 - 2a| 


211 -a| 


If Eqs. (|B.l|) and (|B.2|) the following auxiliary functions have been introduced: 


efHSr) = 


- + a{l-2a) 
^ + p(l - a) 


Hi^\9r) - (1 - a)gr Hl^^^{gr), 


(B.3) 

(B.4) 

(B.5) 


furthermore, following the standard notations, we have that J^l^\gra) = J^a^'^*i9m) and 
= Gp^'^*i9r)- Using Eqs. (B.4) and (B.5) the relation |d+p — |(i_p = 1 can be 


explicitly verified. As in appendix ^ it is useful to investigate the specific limit a —)■ 0, 
ui —)■ 1 and /3 —)■ 1. In this case, defining xi = kri and X 2 = kT 2 the exact form of d±{xi, X 2 ) 
can be written as: 


d+{xi,X2) = 

d+{Xi,X2) = 




IQxi 

^-3ix2 

16 xf X 2 


—2ix2 


+ e^'^‘[8x2 + 4*3-2 - 1][23 i(3i +*)-!] k 


[832 - 4*32 - 1] + e^“‘ [1 - 23i(! + 3i)l 


(B.6) 

(B.7) 


The expressions of Eqs. (B.6) and (B.7) can be expanded in order to derive approximate 
expressions of the transfer function of the relic graviton spectrum across the radiation-matter 
transition. The exact results for the mixing coefficients have been used as a systematic 
cross-check: the power spectra and the energy density of the gravitons in the case of a 
dynamical refractive index are complicated functions of a which must anyway reduce to 
known expressions in the a —)■ 0 limit. This is not only true for the exact expressions but 
also for the corresponding approximated mixing coefficients computed in the limits | kxi \ < 1 
and \kT 2 \ < 1. 
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